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A method of 
enhancing and normalising 
x-ray images, particularly 
mammognuns. by 
correcting tiie image for 
digitizer blur, glare from 
the intensifying screen 
and the anode-heel effect. 
The method also allows 
the calculation of the 
compressed thickness 
of the imaged breast 
and calculation of the 
contribution to the 
mammograms of the 
extra focal radiation. The 
correction of the image for 
glare from the intensifying 
screen allows the detection 
of noise, such as film shot 
noise, in the image, and in 
particular the differentiation 
between such noise and 
microcalcifications. 
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X-ray Image Processing 

Introduction 

5 The present invention relates to improvements in the processing of digitised x-ray 

images (particularly of the breast - termed mammograms or mammographic images), and 
more particularly to the enhancement of such images to assist clinicians in making accurate 
diagnoses based on them. 

X-ray imaging is used as a basis for many medical techniques and, in particular, 

10 mammography continues to be the examination of choice for early detection of breast cancer 
in po5t-menopausaI women and is the basis for national screening programmes. 

Early detection of breast cancer greatly improves mortality rates, perhaps by as much 
as 25%. For this reason, mammographic examinations are nowadays performed on about 25 
million women annually in the EC (of which, about 3 million are in the UK), at a cost of 

15 about 3 Bn US dollars per year. This huge cost and the poor accuracy in diagnosis (8-25% of 
cancers are missed and 70-80% of open surgical biopsies turn out to be benign), have led to 
increased interest in applying computer-aided techniques. Funher, pressure for a reduction in 
the screening interval and fcr routine adoption of nvo-view screening (cranio-caudal and 45 
degree medio-lateral) would entail a large increase in the number of mammograms to be 

20 analysed in the screening programme and this makes the development of reliable and robust 
computer techniques vital. 

As an application of image processing, mammographic images pose a tough challenge 
because they have poor signal-to-noise ratio. This is largely because the images exhibit 
complex textures, and because there is a compromise bietween radiation dose and image 

25 quality. Worse, abnormalities appear as quite subtle, irregular, often non-local differences in 
intensity and the images are ine\ itably ciuaered due to superimposition in the image of 
features separated in the breast. Funher, the background varies greatly between different 
breasts, and there is relatively weak control of the imaging acquisition. 

Unfortunately, while there have been proposals for the application of image 

30 processing to mammography, the vast majority have been of limited scope and incorporate 
only general non-mammography-specific image processing considerations. This involves . 
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great dangers. Image smoothing using such standard techniques may make lesions easier to 
locate, but can remove calcifications and spiculations which can be signs of cancer. Edge 
sharpening using standard techniques may appear to improve an image, but malignant lesions 
typically have fuzzy edges whereas benign ones tend to have sharp edges, so this edge 
5 sharpening process can acmally transform an image of a malignant lesion into one that 
appears to a radiologist to be benign. 

Calcifications present particular interest and problems. Localised cancers in the ducts 
or lobules of the breast are often are associated with secretions that thicken or become 
necrotic. These are called calcifications, or microcalcifications if they are smaller. Clusters of 
10 microcalcifications, which appear as small bright white objects in a mammogram, can be one 
of the earliest signs of breast cancer. Thus the identification of calcifications is a major goal 
of screening programmes, though benign calcirlcations are common (for example they often 
occur in blood vessels), and so the distinction between ductal and vascular 
microcalcifications needs co be made if the number of false positives is to be sufficiently low. 
15 This is not a problem for radiologists interpreting mammograms against the background of 
their knowledge of breast anatomy. It is, however, a challenge for image analysis programs. 
Further, dust and dirt entering the imaging system can create artifacts that mimic the 
appearance of microcalcifications and thus can be a cause of false positives both for 
radiologists and automated imaging systems. 
20 Further problems are caused because the imaging process itself introduces a number 

of variables which affect the image and these will be explained below. 

In the accompanying drawings Figure 1 shows a schematic representation of the 
components of a conventional screen-film mammographic system. When a mammogram is 
performed, a beam of x-ray photons 1 from an x-ray tube 2 powered by a generator 4 is 
25 directed towards a breast 3 compressed between compression plates 17. This beam 3 is 

filtered by filter 5 to remove low energy photons and collimated by collimator 7 to the area 
of interest. The beam has a spectrum of energies that is characteristic of the tube voltage and, 
in particular, the material of the anode 9 but the spectmm is independent of the woman being 
scanned and the view taken. The intensity of the beam exiting the breast is related to the 
30 thickness and type of tissue in the breast. The x-ray photons leaving the breast normally have 
to pass through an anti-scaner grid 1 1 before they reach a phosphorous intensifying screen 
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13. If an x-ray photon is absorbed in the screen 13, light photons are emitted by the phosphor 
and these light photons expose a film 1 5 which is processed to produce a mammogram. The 
exposure to the breast is stopped once an automatic exposure control 19, positioned under a 
section of the breast, has received a set exposure. To generate a digital image the x-ray film 
5 is typically digitised using a laser scanner system or CCD and light box (not illustrated). 

The intensity of radiation incident on the breast in such a system varies spatially for 
several reasons. The most significant is the "anode heel effect". An x-ray tube produces 
x-rays by firing an electron beam at an anode. As the electron beam penetrates the anode the 
electrons are absorbed at varv'ing depths and the x-ray photons that are produced have to 

10 travel through varying thiclcnesses of anode material before leaving the anode. This leads to 
varying attenuation of the emergent x-ray beam thus giving spatial variations in the incident 
x-ray spectrum, this is termed the anode heel effect and is quite substantial. Another source of 
spatial variation is due to the diverging nature of the beam. This means that the further away 
from the source the more spread out the x-ray beam is. However, this effect is small given 

15 that the distance from the source to the breast is large relative to the breast size. In visual 
assessment of mammograms the clinician mostly considers local variations in intensity, so 
the smooth change caused by the anode heel effect is not too troublesome. However, it does 
cause a problem for automated systems. 

Two further effects of the image forming process which affect the image are scatter 

20 and extra-focal radiation. Considering scatter first, the x-ray radiation passes through the 
breast as shown in Figure 2 and what is known as the primary beam passes through in a 
straight line from the anode to the intensifying screen. However, some x-ray photons are . 
scattered in the breast and arrive at the screen from unexpected directions. The anti-scatter 
grid, which typically consists of a series of angled lead strips 21 separated by paper and 

25 aligned with the primary beam, removes many but not all of these scattered photons. Thus 

some scatter will reach the screen and be recorded on the film. Scatter can be estimated using 
the techniques described in "Computing the scatter component of mammographic images", 
by Ralph Highnam, Michael Brady and Basil Shepstone published in EEEE Med. Imaging, 
1994, 13, pp 301-3 13. This allows the calculation of the primary energy imparted to the 

30 screen to be improved by the subtraction of the energy due to scatter. Extra-focal radiation 
refers to radiation which comes around the edge of the collimator 7 as shown by numeral 1 A 
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in Figure 3, and can constitute up to 15% of the total, some of which is scattered and reaches 
the film/screen. Simple techniques for measuring the extra-focal radiation are known, for 
instance from the paper by Highnam, Brady and Shepstone mentioned above. 

The processes of intensification by the uitensifying screen introduces blur or glare 
5 into the recorded image because the absorption of an x-ray photon at a point site 13 A in the 
screen results in the approximately isotropic emission of light as shown in Figure 4 which 
results in blurring of the image recorded on the film. 

Further, the relaticaship between the density of the image on the film and the energy 
imparted to the intensifying screer: is not linear and changes with film processing conditions. 
10 Again, this may not affeci visual assessment which is based on local variations, but would 
affect automated analysis, especially if the non linearity were not explicitly taken into 
account. 

Finally, the process of digitizing the film introduces digitizer blur into the digital 
representation of the image. 
15 It will be appreciated, therefore that the enhancement of x-ray mammograms and 

automated recognition and differentiation of features in mammograms is a very difficult 
problem. 

It has been proposed that the conversion of a digitised mammogram into a particular 
representation, termed the representation, is capable of improving the enhancement and 
20 analysis of such mammogram.s. This was described in "A representation for mammographic 
image processing" by Ralph Kighnam, Michael Brady and Basil Shepstone published in 
Medical Image Analysis; 1996, vol. I, no. 1, pp I -18 and, since the present invention i? 
concerned with improvements to it, will be briefly explained below. 

The intensity of a mammogram at a given pixel (x, y) indicates the amount of 
25 attenuation (absorption and scattering) of x-rays in the pencil of breast tissue vertically above 
(x, y) on the film. 

Ideally, one might hope to be able to produce a quantitative three-dimensional 
representation of the breast with each voxel labelled with a tissue type, such as: glandular, 
fibrous, cancerous, fat, calcium. Given the x-ray attenuation within a voxel it is certainly 
30 possible to classify fat since it has relatively low linear attenuation coefficients. It is also 

possible to classify likely occurrence of calcium, which is practically radio-opaque. However, 
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the remaining breast tissues are those that comprise anatomically significant events in breast 
disease, such as cysts, malignant masses, fibroadenomas, and they are difficult to resolve 
from x-ray attenuation measurements alone. In the h^^, representation these remaining tissues 
are classified as "interesting tissue". Further, there is actually very little calcium so for 
5 practical purposes it can be ignored. 

Unfortunately, a further problem arises because of the projective nature of 
mammographic imaging: the three-dimensional information is lost. In light of this, the only 
information that is available describes the tissue within a cone of the breast, where the cone 
has as its base the area of a pixel and as its apex the x-ray source. After appropriate correction 

10 the x-ray beam within this cone can be considered as a pencil beam. Thus in the /7;„r 

representation (with calcium ignored) there are basically only two tissue classes of fat and 
interesting tissue to consider, and the thicknesses of the interesting tissue (/!,„, cm) and fat 
(hfat cm) which together rr^ust necessarily add up to the total breast thickness H ( i.e. H = hnt 
hfa^ are used as quantitative breasi measurements. 

15 In practice /7,„, is computed from a mammographic image using data related to system 

calibration and image calibration. The x-ray tube output spectrum is assumed to be relatively 
stable but the anode heel effect is corrected for. The mammographic imaging process has 
several parts which might vary from day to day. In order to effect meaningful image analysis 
by computer, it is necessary to know these variations in order to make the images conform to 

20 a standard. To achieve this requires calibration data. The film-screen response, film processor 
and film digitizer are caiibrated by collecting the following data: 

1 . A step wedge fLtm: A film is produced v/ith a stepped wedge made of lucite 
placed along the back of the film and a lucite block placed over the automatic exposure 
control. This film allows us to calibrate the film-screen system and film-processing so that 

25 energy imparted to the iL::2nsifying screen can be related to film density. 

2. A "blank " film: A film, is taken with a short time of exposure with no object 
(breast) present. The exposure has to be short so that the film does not saturate. An exposure 
of 0.04 seconds, at 1 00m.-. and 2SkV for example produces a film that has film densities that 
vary beDA^een 1.8 and 2.6 (despite looking black). This film provides information about the 

30 spatial variations of the incident radiation intensity. 
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3. The digitized image of step wedge film: The film density on each step of the 
wedge is measured so thai: once digitized, the relationship between pixel value in the digital 
image and film density in the corresponding area of the film is known. 

As well as calibrating the system components, data specific to each mammographic 
5 examination is needed. In particular: 

The tube voltage kV); 

The tube current mA); 

The time of exposure (/^s); 

The breast thickness {H cm). 
10 Most of this information is readily available but measuring the breast thickness H is 

currently awkward since the radiographer has to measure it using a ruler; though newer 
machines are incorporating automanc measurement of breast thickness. 

Given a mammographic image, :he thicknesses of interesting and fatty tissue between 
the x-ray source and each pixel can be found by considering the energy imparted to the 
15 intensifying screen at each pixel which is obtained from the pixel values in the image using 
the calibration data. Let E.j,/x, y) be the energy imparted to the screen in the area 
corresponding to the pixel (x, y), L^_,/x, y) contains both scatter and primary components. 
The primary component £/x, y) is determined by subtracting a scatter estimate from the total 
energy imparted as menticned above. 
20 Now for a pixel with h:^^ cm of interesting tissue and h.-at cm of fatty tissue above the 

corresponding area of the intensifying screen, the total attenuation at any energy £ is 
expected to be;: • 

= h i:;<^, V)(ii iaU;^) - pifai{E)) + H fifai{E\ (1) 

25 where the substitution lyjx. y) - H - /h^/x, y) is made. 

In this case, the energy expected to be imparted to the intensifying screen by the 
primary photons is: 

Epix,y-)= (f)iVu.6e,x,y)AM t ;V-''( £-)£5(£:)C(£) x e-^-^^^'^'-'g-'^'^^-'-^^dE , (2) 
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where (j> is the photon fiux for an x-ray :ube voltage of this varies across the image due 
to the anode heel effect; is the pixel area; t, is the time of exposure; Nq^^(E) is the 
relative number of photons at energy £; S(E) is the absorption ratio of the screen to primary 
photons of energy E\ G(E) is the transmission ratio of the grid for primary photons of energy 
5 E\ Ur^c(E) is the linear atrenuacion coefficient of the (typically) lucite (compression plate) at 
energy £; and Z?^/^,, is the - hickness of the compression plate, all of which are known from die 
calibration or image conditions. 

Note that after substiruting Equation (1) into Equation (2) the only unknown is hiJx, 
y). This can be found by equating the primary energy found in the practical case (i.e. 

10 measured from the image; with the theoretical value (i.e. the expected value calculated above) 
and solving the resulting nonlinear equacion to determine /7,„/r. y). 

This process of converting the image into the representation can be visualised as 
converting the original irr:age so that the fat has risen to float on top of the interesting tissue 
surface, then the fat is peeled off leaving the representation hi„f(x, y). Informally, this 

15 representation can be viewed as a surface and clinically significant effects such as masses 
appear as features on this surface, eg. sn:all hills, as seen for example in Figure 5. Note that 
this is fundamentally different from regarding the intensity image as a surface, since the hi„, 
representation is a quantitative measure of anatomical tissue which is distributed through 
venical pencils of the breast. The importance of stems from the fact that it factors out the 

20 imaging parameters particular to the e.Kamination to yield a representation of the intrinsic 
anatomy that is ultimately what is relevant for diagnosis. 

Summary of the invention 

25 The present invention is coucerned with improvements in the calculation of the /7,„^ 

representation, which not only improve che accuracy of the /?y„, representation but also provide 

some very useful results along the vvay. 

A first aspect of the in\'en::on provides a method of correcting a digital representation 

of an x-ray image for degradation caused by a digitiser, the digital representation being a 
30 measurement of the image density for each of a plurality of pixels of the image, the 

measurement being obtained by iliuminating the x-ray image, measuring the attenuation of 
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the light by the x-ray image and calculating from the attenuated light values of the image 
density, 

the method comprising the seeps of: 

calculating from the pixel values in the digital representation the value of the intensity 
5 of said attenuated light corresponding to each pixel; 

applying to the intensity values a modulation transfer function to correct for the 
degradation introduced by the digitiser; and 

converting the corrected intensi:y values back into values representative of the image 
intensity. 

10 A second aspect of the invention provides, in a method of x-ray imaging using an 

intensifying screen to receive x-rays and emit light to be recorded on an x-ray film, a method 
of calculating from the x-ray image density the energy which was imparted to the intensifying 
screen, comprising the steps of: 

calibrating the x-ray film and in:ensifying screen by measuring the response of the x- 
15 ray film and intensifying screen to a plurality of different intensities of received x-rays; 

fitting a theoretical model of the expected response to the measured response, 
and using the fitted theoretical r.iodel to calculate the imparted energy from the x-ray 
image density. 

A third aspect of tlie invention provides, in a method of x-ray imaging in which an 
20 intensifying screen is used to receive x-rays and emit light to be recorded on an x-ray film, 
the image recorded on the film is digitised to give a pixelised digital representation of the 
image density on the film, a method of enhancing the digital representation of the image to 
remove the contribution thereto of glare from the intensifying screen comprising the steps of: 
converting the digital represenintion of image density into a representation of the 
25 energy imparted to the screen; and 

deconvolving the represencaric r. of the energy imparted to the screen using a 
weighting mask defining tlie poinc spread function for the intensifying screen to remove the 
contribution thereto of glare from the intensifying screen. 

A fourth aspect of the inventicr. provides a method of calculating from a 
30 mammogram the compressed thickness of the imaged breast, comprising the step of 

delimiting in the mammc^ram the regiv.n corresponding to the part of the breast which is 
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compressed from the regioa corresponcing to the uncompressed breast edge by detecting the 
smoothness of curves of equal intensity in the mammogram. 

The aspects above may be combined with other known steps to improve the 
production of an A;„, representation of the breast. 

Another aspect of the invention provides a method of detecting microcaicifications in 
a breast from a mammogram of the bre -.st, comprising the steps of processing the 
mammogram to produce an h^^^ represe.'uation according to the above methods, converting the 
value of hi^j for a candidate region into a value representative of the volume of interesting 
tissue and thresholding the volume values to detect as microcaicifications areas of the 
mammogram where the voiume value exceeds a threshold. 

The invention also provides apc -iracus for carrying out the above methods. Further, 
the invention also provider a computer program, which can be provided on a computer- 
readable storage medium, for conrrollirg a computer to carry-out the methods. 

List of drawings 

The present invention will be funher described by way of example with reference to 
the accompanying drawings in which 

Figure 1 is a schematic represer. ration of a conventional screen-film mammographic 

system; 

Figure 2 is a schen^.atic view of a typical anti-scatter grid; 

Figure 3 is a schematic view of a.collimator showing the path of extra-focal radiation; 

Figure 4 is a schematic view of an intensifying screen and film; 

Figure 5 shows an surface of a mammogram; 

Figure 6 shows a film-screen characteristic curve; 

Figure 7 shows an e.xampie of fizting the characteristic curve; 

Figure 8 shows an example of ti.e effects of compensating a mammogram for 

intensifying screen glare; 

Figure 9 shows :wo aim c rasic; profiles with confirmed shot noise in them; 

Figure 10 shows a:: original en^i gy profile and the energies after glare compensation; 

Figure 1 1 shows an exam;::e of loise detection in a phantom; 

Figure 12 shows aa example of .loise detection on a real mammogram; 
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Figure 13 shoes a funher example of noise detection on a real mammogram; 
Figure 14 illustrates the detection of low contrast calcifications; 
Figure 15 illustrates an image sample containing shot noise; 
Figure 16 schematically illustrates a cross-section across a compressed breast; 
5 Figure 17 is a schematic of a cranio-caudal mammogram; 

Figure 18 shows examples of the breast edge; 
Figure 19 shows examples of the breast edge with varying H; 
Figures 20 and 21 illustra:e roughness measures used in the estimation of breast 
thickness; 

10 Figure 22 shows the coordinate system used for extra-focal compensation; 

Figure 23 illustraces the thickness of breast tissue through which extra-focal radiation 
travels; 

Figure 24 shows two examples of the extra-focal component of the radiation. 
Figure 25 shows me moaulation transfer functions of some film digitizers; 
15 Figure 26 illustraies the modulation transfer functions of intensifying screens for 

different values of the attenuation coefficient for light in phosphor; 

Figure 27 illustrates the modulation transfer functions of intensifying screens for 
different values of the atrenuaticn coefficient for x-rays in phosphor; 

Figure 28 are graphs sho wing the point spread functions for two different screen 
20 thicknesses; 

Figure 29 illustrates the rr.odulation transfer function of two intensifying screens of 
different thicknesses; 

Figure 30 illustrates the point spread function of two intensifying screens for different 
gaps between the x-ray film and the screen; 
25 Figure 3 1 illustraies ihe modulation transfer functions of two intensifying screens for 

different gaps between the tlim and screen; and 

Figure 32 illustraies the geometry for calculation of the point spread functions. 

Description of preferred eiubcdiments 

30 

An embodiment of che invention will now be described by going through the steps 
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necessary to calculate the h-^^, representation. These can be summarized as follows:- 

(1 ) Convert pixel value P(x, y) to film density D(x, y) using the digitiser 
calibration data; 

(2) Remove digitiser blur; 

5 (3 ) Convert film density D(x, y) to energy imparted to intensifying screen £*'"Vr^. 

y) using film-screen calibration data; 

(4) Compensate E"'%Jx, y) for intensifying screen glare; 

(5) Compensate Ej^^'^^Jx, y) for the anode-heel effect and diverging x-ray beam; 

(6) Estimate the tnickness of the compressed breast; 
10 (7) Estimate the scattered radiation i^'^'^/r. y)\ 

(8) Estimate the extra-focal radiation y) components; 

(9) Compute the measured primary energy 

£:7(^, v)= ^S^(^,y)- ET\x^y)-ETix^y) 

(10) Compare this measured primar\' energy with the theoretical primary energy 
15 calculated by equation (2) for different values of ^U^. y). Equating the measured and 

theoretical values using a numerical approach solves for hi^/x, yj. 



Digitisation 

20 Considering steps (1 j and (2), mammographic images are digitised by illuminating the 

mammogram with light and measuring the amount of light transmitted through it. Various 
arrangements are used such as laser-scanning devices or devices which use a light box to 
illuminate the mammogram and a CCD camera to record the transmitted light. Digitisers are 
designed to record and outpu: the pixel values related to the film (image) density rather than 

25 the transmitted light itself because the film density is independent of the illumination level. 

Step (1) of the method aoove is relatively straightforward because modem high- 
quality laser scanning devices, for instance, have a known linear relationship between the 
film density and the pixel value: 

P(x,y) = mD(x,y) -r c , where m and c are constants. 

30 Thus D(x. y) can easily be found from the pixel value. 
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A first aspect of the invention concerns step (2), the removal of digitizer blur. The 
modulation transfer functions which show the degree of blur for various digitizers are known 
or can be measured. For instance, the modulation transfer functions for a scanning 
microdensitometer, CCD camera and a laser scanner are disclosed in "Digital Mammography 
5 - the comparative evaluation of film digitizers" by D. H. Davies, British J. Radiology, 
66:930-933, 1993 and are shown in Figure 25. With the present invention, however, these 
functions are applied not to the film density value, but to the transmitted light 7/ which is 
related to the film density in a non-linear way as follows:- 

. Ti(x,y)= Ii{x,y)x 10'^^"^^^ 

10 where // is the light illuminating the film in the digitizing process and T/ is the light emerging 
from the film. This transmitted light is calculated and the modulation transfer functions for 
the digitizer are applied to it to remove the digitizer biur. Then the transmitted light is 
reconverted into the film density (now without the digitizer blur) by dividing by // and taking 
the log to base 10. 

15 

Film screen response 

The next step, step (3) is to calculate from the deblurred film density D(x, y) the 
energy imparted to the intensifying screen £?"'"^p^/x, y) by using film-screen calibration data. 

20 The film-screen response to energy imparted to the intensifying screen is given by a 

characteristic curve like that in Figure 6. Figure 6 shows a film-screen characteristic curve 
found by exposing a lucite step wedge. Lucite is chosen because of its similar absorption and 
scattering properties to a breast consisting of half fat, half interesting tissue. A small 
correction for scattered radiation can be made, and the logarithm of the calculated relative 

25 energy impaned is plotted against film density, which is measured with a densitometer. 

With the invention a serpentine characteristic curve is fitted to the lucite step wedge 
data, and this fitted relationship is used to calculate the imparted energy from the film 
density. In more detail, the serpentine curve can be expressed as: 

x'y + a'y- b^x - 0 

30 rearranged this is: 



wo 00/52641 



PCT/GBOO/00617 



-13- 



This has the properties: 

Intercept at (0, 0) 



Extrema at (±a,± ) 



a 

Inflection 2X x = Q,x = ± —f= 

v3 



-b' b' 

Extent: -mf < x < -f-inf. — " < V < T— 

2a " 2a 



Symmetry (0,0) - rotational 
Asympcote y=0 



Gradieni 



dx {x'^a") 



dy 

Gradien: at origin = — r 
dx a' 

Thus a,b and also the origin to use need to be detemiined. The origin has to be the 
place of symmetry so there are four unknown parameters: a. b, x^^, J'origm- The equation of 
the curve becomes: 

b (-^ "* X origin^ 



(jC ~ Xorigin^ T a 



Where x = log £ (where E is the energy impaned to the intensifying screen) and y = D 
(where D is the film density). The origin is initially guessed to be in the centre of the 
film-density range and :he fiim density range is measured so as to determine initial estimates 
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of a,fa. An optimizing routine is used to find the optimal values of the parameters based upon 
the step wedge and blank film calibration data collected from the mammography unit. 
Inverting gives: 

JC — JCorigin T — ^— ^— — — — — 

2^' ~ y origin 

5 Figure 7 shows an example of the curve fitting. In Figure 7 the horizontal axis is die 

logarithm of the energy imparled to the intensifying screen. The vertical axis plots the film 
densities. The line marked step wedge data is the mosc accurate and is based upon 
measurements with a piece-wise linear fit between the points. The short dashed line is the 
serpentine fit to that data showing a good fit. BCU is data collected from the Oxford Breast 
10 Care Unit and the longer dashed line then is the serpentine fit to the data collected at the BCU 
which gives a good fit in the low film densities, but not such a good fit at high film densities. 
The fit can be im*proved to give the curve with long and short dashes by noting that although 
the film-screen curve is not synimetric in its entirety, the highest film densities are not used, 
so D,„ax can be cunailed so thai symmetry is obtained. 

15 

Glare removal 

Having calculated the energy impaned to the screen, the next aspect of the invention 
concerns the removal from this energy value of the contribution to it by intensifier screen 

20 glare in step (4) which is directly related to light exposure to the film. In summary this is • 
done using a point spread function calculated by assuming that the x-ray photons are 
absorbed equally across each part of the intensifying screen corresponding to one pixel and at 
different depths within the screen. Knowing the pixel size, the solid angles between each 
potential absorption sit3 and the neighbouring pixels can be computed. These solid angles 

25 define the proportion of light photons that the neighbouring pixels receive. The solid angles 
are weighted by distance from the absorption site to the pixel (representing light absorption 
by the screen) and by tl:e actuiil x-ray energy reaching that site. Using these values glare can 
be estimated and removed froi^- an image. The mathematics of this will now be explained in 
detail. 
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To derive the point-spread function, the screen is considered to be divided into layers 
and then each layer to be sub-divided into sub-pixel size units. Each of these units is 
considered to be a potendal si:e of absorption of an x-ray. Let tp cm be the thickness of the 
intensifying screen. The screen is considered to be in n layers so that each layer has a 
5 thickness of cA^= t/n cm. The layer is subdivided into pixels (whatever resolution is being 
used) and then each pixel is split into 100 smaller elements. It is assumed that the x-ray film 
lies directly on top of the intensifying screen separated only by a screen overcoat so that there 
is a gap of /ocm between the x-ray film and the intensifying screen. For each layer a 
weighting mask w.(x, y) is computed (effectively the point spread function for that layer) 
10 which gives the percentage of light photons emitted at (x^, yj, depth r, that reach the film 
corresponding to the spadal position (x, y). The proponion of photons from (x^ yj reaching 
(x, y) is related to the solid angie d9 from fx^ yJ, depth z to fx, yJ, Symmetry round the 
azimuthal angle is assumed and the ID case considered (see Figure 32): 



15 



dO = cos 



'I 



2\ 



Where b, c are found from simple geometry and Xj and are the distances to either 
side of the target pixel (x, y) from the central pixel (x^ y^J, 

Some of the light photons that are emitted from (x^ y^ z) are absorbed by the 
20 phosphor so that cl9 is weighted using Beer's law and the relative glare becomes: 



relative glare = dOe 

Where fJ^fhosphcr an average linear attenuation value. Dividing the relative glare by 

the total relative glare gives a v/eighting mask for each layer. The x-ray energy being 
25 imparted to each layer is now incorporated. The energy into each layer is: 
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where E,^"*^ is the energy imparted and E"" is the incident energy. So the relative glare now 
becomes: 

r* pnospnor 

relative glare =c/^6 C'C^^oX) 

These values are combined for each layer and for each sub-pixel to get the full 
5 weighting mask w(x, y) and the results are scaled so that I^^yy>i<x. y)-\S) and rotated around 
the central pixel to give a tv^o dimensional mask. 

The parameters ju^p'Z::hor ^ M^^hor included in the model more for illustrative 

purposes than practical value since it turns out that they have little effect on the PSF, which 
is the crucial factor. Figures 26 and 27 show the modulation transfer functions (MTFs) of 

10 screens in which the varying parameters are u'^'iZ^j.^r MpfJsphor • The modulation 

transfer function is simply the Fourier transform of the point spread function. Figure 26 
shows the MTF of intensifying screens using the model with 100 micron screen thickness 
and 0 gap between the film and the screen. There are three curves representing three 

different values of f^^p^Jsphor ■ -"^^ -^^^ value rises the more direct path to the film 

15 becomes more heavily favoured and so the MTF improves. Figure 27 shows the MTF of 

intensifying screens using the model with 100 micron screen thickness and 0 gap between the 

film and the screen. There are three curves representing three different values of f^^pHo^hor • 

As the value rises more x-rays get absorbed nearer the film and so the MTF improves'. 

The two most imponant parameters in the model are the screen thickness tp and the 

20 gap between the film and the screen Figure 28 shows the PSFs for two different screen 
thicknesses. The graphs show the percentage of light photons emitted from the central pixel 
(where the peak is ) and travelling to the surrounding pixels. The graphs are profiles across 
the 2D masks for when the spatial resolution is 50 microns. On the left, the screen is 100 
microns thick, whilst on the right its 200 microns thick. The thicker the screen the greater the 

25 absorption of the x-ray photons at depths further away from the film and thus the more spread 
out the light emitted by the screen. The MTFs (for zero gap) are shown in Figure 29. Figures 
30 and 3 I show the PSFs and MTFs of when the gap between the tllm and the screen is 
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varied. Figure 30 sho vs the PSF of rvvo intensifying screens and using the model with 100 
micron screen thickness and C and 100 micron gap between the film and the screen as 
marked. Figure 3 1 shows the MTF of two intensifying screens using the model with 100 
micron screen thickness and 0 and 100 micron gap between the film and the screen as 
5 marked. Clearly, these parameters have a large effect on die screen performance and they are 
an important part of the calibration data. For example the intensifying screen Fiji UM 
Mammo Fine enclosed within a Fuji EC-MA cassette has a thickness (/p) of approximately 
100 microns. The gap between film and screen is much more difficult to estimate and 
information from manufactur:^rs is sparse. Ffowever, there is probably at least a 1 micron 
10 protective coating on the film emulsion and a funher protective coating of 5-15 microns for 
the screen. 

To properly compensate for glare it is necessary to know exactly where the edges of 
the mammographic film are since no glare comes from outside die film and those areas 
should be treated as zero in the convolution. Also, the energies imparted outside the breast 
15 area but on the film saairate the film so that to model the true effect of glare the energies on 
those regions have to be set tc rhe expected incident energy as computed using the known 
time-of-exposure. The weighring mask w defines the point spread function for the 
intensifying screen so ihat the energy impaned image attained is simply the result of the 
energy imparted without glare convolved with w: 

20 £'giarc ^ £...o glare 

This can be solved by deconvolution in the Fourier domain. Figure 8 shows an 
example of the effects of compensating for the glare. 

Noise detection 

25 

Although the glare ha.- been removed to enhance the calculation of in steps (5) to 
(9) above, this step is itself imponant in that it allows the detection of microcalcifications and 
in particular allows a distincii jn to be made between microcalcifications and film-screen 
"shot" noise which can look c:nfusingly similar to an image analysis system. 
30 Film-screen "shot" noise can arise from dust and dirt on the intensifying screen or 

from deficiencies in the film. A major difficulty in detecting microcalcifications is that this 
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noise tends to appear with simiiar characteristics to the calcification: small, low film density 
(bright) and high frequency. Consequently, automated detection of microcalcifications as 
localised bright spots rends tc generate many false positives. Although some of these can be 
eliminated by using the clustering property of real microcalcifications it would be preferable 
5 to be able to eliminate them individually. The present invention uses the absence of blur from 
bright spots to mark them as noise. The absence of blur indicates they were introduced into 
the imaging chain after the glare from the intensifying screen. 

In the above method c: glare removal the mammogram is first transfomied into an 
image which represenis the energy imparted to the intensifying screen. In this representation 

10 calcifications and noise have very' low values due to their apparent "high x-ray attenuation". 
When glare is removed from, the image, pixels with low energy values due to noise become 
negative indicating that the original energy value was not feasible. Figure 9 shows two film 
density profiles with confirmed shot noise in them. It shows two profiles across 100 micron 
images containing microcalcifications. The pixel value on the y-axis is linearly related to film 

15 density, with high pixel value meaning low film density and thus low energy imparted to the 
intensifying screen. In the lef: profile the image restoration algorithm for glare picks out r=50 
as being film-screen noise, and in :he right prorlle the point at x=206. The peaks at>:=190 
and x=250 in the right profile are real calcifications. Figure 10 shows an original energy 
profile and the energies after ^lare compensation. It shows uvo profiles across a 100 micron 

20 energy image. The left profile shows the energies imparted before glare compensation. There 
is a low value at .\:=62 which is due to noise (low energy means it appears white). The right 
profile shows the energies af.er glare compensation, the energy at the noise point is now 
negative, i.e, it has been introduced to the image after the major blurring stage. 

Figure 1 1 shows an example of a physical phantom (TOR-MAX, University of Leeds, 

25 Department of Medical Physics, X-ray Test Object) and the noise that the detection scheme 
finds. The image on the left is of a physical mammographic phantom. The phantom has 
simulated microcalcirlcations (pointed at by the arrow), round cylinders of var>'ing densities 
and, at the far right, patches of texrure. Before the mammogram was performed the 
film-screen combination was opened up and dust and dirt distributed randomly - these 

30 show-up as being bright whr.e on the left image. The results of the noise-detection algorithm 
are shov/n on the right. All ::-e obvious noise pixels have been detected as well as some far 
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more subtle ones, none of the simulated calcifications have been marked as noise. 

Extensive testing of the filni-screen shot noise detection algorithm has been carried 
out on real mammograms. On a set of 20 sections of mammograms many of which contained 
microcalcifications an experienced radiologist marked the points corresponding to 
5 fiim-screen shot noise. She marked 1 56 points of which the above algorithm detected 150. 
The algorithm also detected 6 more points which the radiologist was unable to state 
categorically whether they were noise or real microcalcifications. No definite 
microcalcifications were marked as noise. Figures 12 and 13 show examples of mammograms 
and the noise detecred. Referring lo Figure 12 on the left is a section of a real mammogram 

10 digitised to 100 microns per pixei. There is a cluster of noise points as confirmed by an 

experienced radiologist and by the algorithm as seen on the right which shows in white the 
noise detected. There is a quire obvious hair in the top right corner and some subtle genuine 
microcalcifications in the bottom left comer which have not been marked as noise. Referring 
to Figure 13 on the left is a section of a real mammogram digitised to 50 microns per pixel. 

15 There is a whole cluster of noise points as contlrmed by an experienced radiologists and by 
the algorithm as seen on the right which shows in vvhite the noise detected. 

Although this noise deieciicn scheme works extremely well, there are two 
considerations which should be noted: 

(i) At points where the film is saturated, that is where it has very high film 

20 density, or at very low film density such as in those areas beneath lead markers the estimated 
energy imparted is inaccurate and that can affect the noise detection. However, this happens 
only at the very edges of the breast image well away from any likely calcifications. An 
example of this is in Figure 1 1 where the very low film densities at the top of the image are 
marked as noise. 

25 (ii) Noise points affec: the noise detection at points near to them. The glare 

removal requires the energies in a local neighbourhood. If one of the energies is artificially 
low due to noise then the glare removal might be incorrect and might cause other noise not to 
be detected. This can be seen in Figure 1 1 where there is a hair which clearly has the end 
points marked but not the central ones. This can necessitate a two-pass glare removal 

30 algorithm or require that any pixels nominated as calcifications within a certain radius of a 
noise point are also considered noise. 
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This method for detecting microcalcifications can be further improved. It may be 
recalled that the fundamental assumption underlying the generation of the representation 
is that the breast consists entirely of fat and "interesting tissue". Since calcifications have an 
x-ray attenuation coefficient which is about 26 times higher than those tissue types, the 

5 attenuation of an x-ray beam through a microcalcificaiion, perhaps of diameter 0.5 mm, is 
comparable to that through 1 .3 mm of interesting tissue. For this reason, the value 
computed for pixels which in fact correspond to calcifications is expected to far exceed those 
which correspond to non-calcifications. 

Based on this the improved method uses the "interesting tissue volume", denoted by 

10 v,„„ which represents the total amount of interesting tissue present in a mass of breast tissue. 
Vf„, can be computed from the lh„, values over the region of interest on the image. Now 
consider the v,„, value of a small volume of breast tissue, 5, whose actual volume is v^. If no 
calcification is present in B, v,.„, should be bound above by v,e,. However, if 5 is a 
calcification, the computed v, „, would exceed v^, owing to the violation of the fundamental 

15 assumption of the ih,„, model. Thresholding the V:„, to v^, ratio, enables differentiation of 
calcifications from other breast tissue. This ratio is subject neither to varying imaging 
conditions nor to diftsrent tissue backgrounds on which a calcification is projected. This is 
contrary- to image conirast, which is what most other calcification detection algorithms use. 

Candidate regions of the mammogram can be segmented e.g. those that satisfy a weak 

20 contrast constraint and which are not too large. The calcification, of course, has other breast 
tissue above and below it which because of the projective nature of the mammogram 
contributes to the. image. This contribution is removed by estimating the background. /i,„, 
(from the /?,„, values surrounding the candidate region) and subtracting it. v,„, can then be 
calculated. The estimation of v^., is more difficult, and to do this the heuristic assumption is 

25 made that the candidate microcalcification has an elliptical cross section, so that its volume 
can be estimated from its projectioii in the image. 

In an initial experimental study a total of 20 image samples taken from 7 different 
mammograms were used. The images are digitised to a resolution of 50|j.m per pixel. There 
are altogether 27 microcalcincaiions in the 20 samples, each of which has at least 1 

30 microcalcification. 100% true ;;0sitive rate is obtained along with 0 false positives per 
imaae when the v,,J-'' threshold ii set to 3. The algorithm detects 4 other regions which the 
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radiologists are unable to make conclusive remarks on whether they correspond to real 
calcifications or non-calcifications. 

The results were also compared with those obtained by simply thresholding the grey 
level contrast. An ROC analysis shows that this method achieves both higher sensitivity and 
5 better specificity, albeit on a small sample. Figure 14 shows an example which compares a 
typical contrast-based detector attempting to detect a low-contrast microcalcification with the 
improved method above: 

(a) the original image with 2 microcalcifications; 

(b) detection result using the volume ratio algorithm: both calcifications are 
10 detected successfully; 

(c) detection result using grey level contrast: the more obscure calcification on the 
right is missed. 

Thus the glare removal process of the invention enables the elimination of many false 
positives that correspond to noise. Figure 15 shows an image sample containing image shot 
15 noise: 

(a) the original image vvith 1 microcalcification and a number of image shot noise 
points scattered over the image, of which the one lying horizontally near the bottom of the 
image is very prominent; 

(b) detection result using the volume ratio algorithm: the region corresponding 
20 to shot noise is eliminated and oniy the real calcification is marked, 

(c) detection result using grey level contrast: the image shot noise is wrongly 
marked as a calcification. 

Algorithms based on contrast measures have difficulty rejecting such false positives. 
Thus despite the h,^^, representation explicitly calcifications it is useful for detecting them and 
25 for differentiating false-positives due to shot noise. 

Anode heel effect 



Remming to improvement in the calculation of /?,-„, , the next step is to correct for the 
30 anode heel effect in step (5). 

The primar;/ component of incident energy at (x, y) is directly proportional to the 
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number of photons incident to the volume of tissue projected onto that pixel: 

^7(^, V) = j.)^^^^7^^/^,y) , (3) 

where the last term is used to denote that part of Equation 2 which is independent of the total 
number of photons (nd stands for not-dependent). It is assumed that the x-ray energy 
5 spectrum stays the same but that :he total number of photons ^ (V„ x, y) changes with (x, y) 
due to the anode heel effect and diverging beam. The scatter component at the pixel (x, y), 
mostly comes from the x-ray photons that are entering the breast tissue in the surrounding 
neighbourhood. This neighbourhood is small enough to allow the anode heel effect over it to 
be ignored so that the scatter component is also directly proportional to the incident radiation 
10 at^x,;^;: 

The total energy imparled is the sum of the primary and scatter components so that 
using Equations 3 and 4 gives: 

15 The incident photon tlux is greatest underneath the anode, let the position on the film 

at this point be (x^ yj. The aim is to change E^^fx, y) to be as if from that incident photon 
flux. 



SZr.c.M,y) = -777r7^E""''{x,y), (6) 



The ratio of the two photcn fluxes needs to be computed. To determine the anode 
20 heel effect for a specific system and to thus compute the ratio an x-ray exposure is performed 
with no object present which gives an apparently "blank film". For example, nine points on 
the film can be sampled. The energy impaned to the screen comes mostly from the primary 
radiation since there is no scattering material: 

25 ^Z.(-..X:) = ?(K.^.,yJ(sAAE7Jx^,yJ+ Er,^ix„,y.)) 
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Noting that since the pans of the signal not dependent on the number of photons are 
equivalent (there is no object) gives: 

<l>{V„x,y) ~ EZ,ix,y) 
Substituting this into Equation 6 allows compensation for the anode heel effect. 

5 

Estimation of compressed bre:-st thickness 

.-^n important step in this method is the estimation of the compressed breast thickness 
H. This is needed for an accurate calculation of but is also useful in checking the 
10 accumulated radiation dose per unit volume, and in locating the nipple which is important for 
diagnosis. 

Many recent mammography systems have built-in analogue or digital thickness 
meters but their accuracy and precision are currently wanting. In clinical practice, most 
existing systems do not have such indicators. 

15 This aspect of the invention provides a robust and accurate method for estimating the 

compressed breast thickness from a mammogram using image processing and modelling 
techniques above. The estimation is based upon the e.xistence of the "breast edge", a 
fatty area around each breast where the breast thickness steadily reduces to zero. Determining 
that area using image processing provides enough data to estimate the breast thickness when 

20 calibration data such as the tube voltage and exposure time are known. The technique can 
also be applied to mammograms caken previously. 

During mammography, the breast is compressed between two supposedly parallel flat 
compression plates. The com.pression causes the breast to spread out, so that over most of the 
plate the breast is of equal thickness. However, towards the edge of the breast the breast 

25 bulges like a balloon and there is a not a straight venical edge; Figure 16 shows this 
schematically. This bulge is v.'hac we term the breast edge. 

The breast tissue is enveioped in two layers of fibrous tissue, the deep layer overlying 
the muscle, and the very thin superrlcial layer beneath the skin. The superficial layer is 
separated from the skin by 0.5 to 2.5cm of subcutaneous fat or areolar tissue. Joining the 
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layers to the skin are fine fibrous ligaments (Coopers ligaments). From this, it is reasonable 
to assume that what is termed here the breast edge consists entirely of fat, though in localized 
regions, particularly near the nipple and near the ligaments, this will not be strictly true. 
Figure 17 is a schematic of a cranio-caudal mammogram (not to scale). The left picture 
5 shows the breast, while the right picture is the mammogram. The shaded areas are meant to 
represent volumes of interesting tissue. The arc denoted by the letter E lies right on the edge 
of the breast as seen from the x-ray source, that is, an x-ray from the source tangentially 
touches the breast before reaching E, There is very little attenuation of the x-ray beam 
anywhere along E. Along the arc dep.oced by the letter D, there is a greater thickness of tissue 
10 than at E, but still relatively litiie attenuation, since this is the breast edge and most of the 
breast tissue is fat. Along the arc C, there is H cm of fat so that the x-ray attenuation is quite 
large but it is still uniform along the arc. This changes for arcs A and B, where there is a 
heterogeneous mixture of tissues. The x-ray attenuation along these arcs can vary from being 
due to H cm of fat to almost H cm of interesting tissue. 

15 This method of determining breast thickness H is based-upon delimiting the projected 

breast edge from the interior of the breast, i.e to determine the arcs C and E in Figure 17. 
The breast can be segmenied from the background to determine the smooth arc E on the basis 
of the film densities. Then ail the pixels which are slightly brighter than those on E are 
marked to get another smooth cur^-^e, since this is still in the homogeneous breast edge. This 

20 may be arc D. Eventually, arc C is reached, which is still smooth, since it is just within the 
homogeneous breast edge region. If pixels are continued to be marked which are slightly less 
dark, it starts to mark pixels within the interior of the breast. These could be anywhere 
and numerous since the interior is a heterogeneous mix of tissues. Thus the arc that delimits 
the interior of the breast (arc C) is -he extremum of the smooth iso-intensity curves. From this 

25 model, it is also clear that the breast - issue between arc C and the x-ray source comprises H 
cm of fat. Thus, if C can be found, and given a good model of the mammographic imaging 
process, the pixel values along C can be used to estimate the amount of fat which must have 
been between the pixels aiong C and the x-ray source thus giving the breast thickness H. 
There are two possible approaches to finding the curve which delimits the breast 

30 interior and then estimating the breast thickness. One is to find the smooth curve by using the 
original pixels values vvhich are Lfoduced fay the digitiser. The imaging process is then 
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simulated using different breast thicknesses until the predicted pixel value matches that which 
is found on the smooth curve. The second approach is to take a rough initial estimate of the 
breast thickness and then to generate the representation. Analysis of the representation 

enables determination of hovv :he estimate of the breast thickness should be changed. When 
5 the /7,„, values are generated using an accurate breast thickness, it is found that within the 
projected breast edge there is so iittie attenuation that there cannot even be H cm of fat within 
that region. For these regions is set to zero and then it is determined what thickness of fat 
alone would give the observed ar.enuation. FigurelS shows some examples of the breast 
edges found by marking those pixel.; that have /7,„,=0. The bright white areas represent the 

10 so-called breast edge. This is where ihe breast stans to decrease from its constant thickness H 
and the attenuation is no longer eriough even to be of H cm of fat. The inner edge of the 
projected breast edge region is maihemadcally quite smooth. Compensation for compression 
place slant (typically up to 0.5cm from the chest-wall out to the nipple) is performed by 
effectively adding a wedge of fat to the breast by adjusting the mammographic image using 

15 monoenergetic modelling assume ticns. 

The method of determiaiiig H is to start with an underestimate of the breast thickness 
and then to compute the reprejeatation. Those pixels that have hi^^-O are marked 

and a measure of how rough the cur/e that those pixels represent is computed (see below). 
Initially, the Z?,-,,; values wili be tco high and there will be no pixels with hi„i-0. As H 

20 increases, more and more pixels are detected v/ith /2;„.=0 and these pixels will represent a 
smooth curve. This continues un::; :he internal breast region is reached at which point the 
roughness measure rises dramaticr.iiy to indicate a rougher curve. 

Performing the computation this way x-ather than just using the pixel values in the 
original image and looking for isc -intensity curves allows more accuracy since H can be used 

25 to predict scattered and extra-focal radiation. It also enables checking that the algorithm is 
working properly. For example strirting with many points that have > ^ then the value 
ofH is far too low and shoaid be ir.creased it dramatically rather than in small steps. 
Another example is that if more en a: 20% of the breast is found as being within the projected 
breast edge, then H has incre:isec tco far. Figure 19 shows real examples of the breast 

30 edge as H varies. This shows imr.2:^:s with luminance proportional to hf„f except for the bright 
white breast edge whi is wh ^rf: . :. .-(). The breast thicknesses tried are marked above each 
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image. The firsr value, 3.4cni, is too low - there is no breast edge; the second value, 5.4cm 
is just about right - there is a subr.aatiai breast edge and a smooth internal edge; the last two 
estimates are far too high - the projected breast edge has become ragged. 

As for the computation of roughness, as breast thickness increases a large jump in the 
5 measure is expected just after the correct value. In Figure 20, the curves marked A-E would 
ideally have the roughness measures as shown on the right. 

There are several possibil:::es for such a roughness measure but the approach in the 
invention is to use an estimate of .i fractal characteristic of the curve. A fractal curve has two 
parameters: fractal dimension anc 3 -dimension, and for estimating image textures the latter 

10 is more reliable, stable imd gives ::e::er discrimination. A technique to estimate the 

D-dimension: "the covering blanriet" which is based upon morphological operations is used. 
To compute the measure, two fui-r.er curves are created from the curve: one from 
opening the /7/„,=0 cur\^e and one rVcm closing it (see Figure 21). An opening operation 
creates a smaller curve which is s:ncother than the original, whilst a closing operation creates 

15 a larger curve which is smoo:her .r.an the original. The area contained between the two curves 
is a measure of roughness: the larger the area the rougher the curve. A 5 pixel by 5 pixel disk 
is used as our morphological operritor and the thickness estimation is performed on 300 
micron resolution images. 

Although the area bervvee:-. tl\e curves is related to roughness, it is also directly related 

20 to the total length of the cup. e, so is normalised and made dimensionless by dividing by 
the length of the curve. Taking a /: t.iese considerations into account the roughness is the 
number of pixels enclosed by closing and opening of hi„,-0 divided by the number of pixels 
on the hi„r=0 curve. When this n^easure increases above a threshold of 1.3 increase ofH is 
stopped. 

25 In this estimate breas: ihickness, the smooth curve which delimits the interior 

breast region is detecre::. The prr>;erce of Cooper's ligaments in the image, which is quite 
unusual, would affect cniy a sec:: :.n of the curx-e and only slightly the overall roughness 
measure. The Increasec density c:" ::3sue around the nipple just means that the curve is not a 
semi-circle but a more irregular £:-.:i;:e that is still smooth. Indeed, the fact that the curve is 

30 less regular can be usee lo decec-. position of the nipple. 

There are sevenii \va\ i :;: :o:ming an initial (under)estimate of breast thickness from 
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the calibration data and image. One vvay exploits the fact that near the chest wall there is low 
scatter and the breast tends to be rar.y so one can assume /7,„,=0. hf^^- H, From the energy 
imparted H can be estin:ated. A:ioi:: .er way is to assume that near the breast edge there is 
pure fat and some nominal, high sea. ler-to-primary ratio. Another way is to estimate an 
5 initial value for the breast thickiiess using a film density that is known to be outside the breast 
and the calibration data; this always gives an underestimate of the actual breast thickness. 

Bounds on the breast thicir.ess can also be estimated to check that the method is not 
trying infeasible breast -hicknesois. The lower bound on H is related to the minimum 
attenuation apparent wi:hin the ::rea..r image. To achieve such low attenuation requires a 

10 certain minimum thicia':ess of breas . tissue. The minimum possible H occurs if the breast 
tissue has very low fat conteni: sc ::iat = H (i.e. highly attenuating). An upper bound on H 
can be determined in exactly :he ?.::rf;e way except using the maximum attenuation and 
considering the breasi to be near/- a'i fat. 

As indicated in :he inircv-ici: on, the scattered radiation can be estimated using a 

15 published technique, achieving v.tz- 7), and then the extra-focal radiation must be removed in 
step (S). 

In this embodiment the evtra-focal radiation is found by using the "inner projected 
breast edge" found in the estirr:a::cr: of// above. 

The extra-focal radiatio:. assumed to be travelling in random directions much the 
20 same as light from the sun on a .: iy day and is arriving at a point on the intensifying 

screen from all directions equaii y. ,f ihe thickness and the composition of the breast along 
some curve C is known, the expected primary and scatter components along C can b.e 
estimated. The extra-focal component along C can be found by subtracting them from the 
known total energy imparted: 

this value is then extrapolated i- ' t:' :.ie entire image. 

As mentioned ajove :he rur e C used is the "inner projected breast edge". This is the 
breast tissue around the edge c: '.le :oir;pressed breast which is fat and which is where the 
breast starts to curve. Figure 3 s .ov. this piciorially. In Figure 3 it can be seen that most 
30 x-rav radiation comes from the. jca spot. However, some extra-focal radiation is scattered so 
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that it reaches the film-screen and con-j-ibutes to the final image. Since the breast edge is 
curved, photons IB arriving at low arigles have easy access to the film-screen and thus have 
a significant impact on the final ir.iag^:. In the following the breast edge is defined to be that 
volume of the breast which is far :.nd 3 v/here the breast thickness descends from H towards 
5 zero. On a mammogram this area can oe depicted by an inner and outer curve 3A and 3B as 
shown in the top-right of figure 3 and ehese can be calculated using the techniques described 
above for the calculation of H. 

The inner curve is used as cur. e C and the breast edge is modelled as consisting 
entirely of fat and being senii-circular. Using this assumption, the thickness of tissue 

10 berween any point 

and where the extra-focal is perceived to be coming from can be computed and then this is 
used to adjust the extra-focal estimate accordingly. 

Figure 22 shows the coordinate system used. (X V, Z) are used as the 3D coordinates 
of a point, with the X axis along the chest wall, the Z axis perpendicular to the chest wall, and 

15 Y the vertical axis rovvards the »;ay lOurce. fx, z) is a point in the 2D image plane. When no 
breast is present the fiim-screen comb ination at (x, z) is assumed to receive the same 
extra-focal radiation from all directions ( 6, p). p is termed the azimuthal angle, 0 the 
elevation angle and E'^^^xira is trie ''extra-focal constant", that is, the energy incident on (x, z) 
from any direction. 

20 The energy imparted to th - scr een due to extra-focal radiation can be assumed to be 

constant across the image if there vver- no breast present. With the breast present, the 
extra-focal radiation is attenuatec by that quantity of breast tissue present along any particular 
path of travel. 

That quantity of breast tissue 13 estimated using a model of the breast edge and the 
25 constant E^^^cxira that would have beer, imparted had no breast been present is computed. 

At each pixel (::, z), the ei-erg-^ impaned when a grid is used can be approximated by: 

30 where T is the transmission throu'-h t::e anti-scatter grid, and assuming that the screen 
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absorbs all the photons reaching it, regardless of angle. To simplify the analysis a 
monoenergetic case with phoior. energy e is treated. The total extra-focal radiation imparted 
to the intensifying screen at any point when no object is present is given by: 

Er(x,2)= \\E7{x,y,9)dddp = EzA\T{9,p)dedp 

0 0 0 

Let /jh(x, 2, 6, p) be the a:ien.:acion due to the breast along the path to point (x, z) from 
angles : 6, p. 

with s the relevant photon energy. The energy impaned due to extra-focal radiation can be 
estimated using Beers Law, ignoring scatter, for when the breast is present: 

Er(x,z) = EZ. ! \%-'"''''-''-''T{9,p)dedp 

To study the problem analytically, .this equation must be simplified. This is done by 
first assuming symmetry around the azimuthal angle p so that the problem becomes 
essentially one dimensional: 

• Er(x,z)=27cEZ.\e-'''-''rie)d0 

E^%.xira is estimated toge:her with the attenuation along each ray using a suitable 
model of the breast and breast si .ape as explained below. The problem can be reduced further 
since 7t/2<9 <k represents the ai.glei of extra-focal radiation coming from the chest wall 
where there is a full thiclcness o: breas: and body, so that contribution is assumed to be 
minimal. Furthermore, the mos : imponant part of the breast for extra-focal radiation is the 
breast edge and that can be ccns-dered =0 be just fat, so that the equation can be further 
simplified: 
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2 



In this equation hfa/x, r, 9), i-^- thickness of fat from any point along any direction, 
is computed as explained below. Using the coordinate system of Figure 22 the shape of the 
breast edge is considered to be semi-circular With no loss of generality save for assumptions 

of rotational symmetry' about th-^ azimuth, the ID case where X=0 and Zis perpendicular to 
the chest-wall is treated, so that the breast edge can be described by: 



r(z,) = ^±V(^7^-zj), 



where a -H/2 andZo is the Z-ccordinaie relative to the edge of the inner breast edge, i.e. Zq = 
Z - Ze^^e- Figure 23 shows the two cases inside the breast. The case Al shows the ray 
10 travelling through both breast and breast edge: for angles of 0 < 9 <6h where 9h = tan"' H/] Zq 
I the extra-focal rays must travel through a thickness of breast tissue and a thickness of breast 
edge tissue so that: 



Z Z 

/ //^v ^0 , -^intersect 



COS 6 ' cos 6 

and Z.:„icrsect is the intersection of the ray and the breast edge which occurs at: 



m(a - c) i J(m'a' - c*" -r lac) 
15 7 = — ^^-^ — 

^intereect 1 . .-,^2 

where m = tan 9 and c - Zq tan 9, For the points which satisfy case A2, some of the 
extra-focal rays will come through tissue which is within the uniform breast thickness and 
some through the breast edge area. These rays are from Oui to tt/Z and the thickness of the 
breast is simply: hf^t(S) = H/'siji 9. Outside the /7y„,=0 line, there are angles that have a free 
20 line-of-sight to the pixel (x, z), ^;ee case Bl in Figure 23. In fact, these angles are from 0 to 
9in„it where ^7„„-, is the angie of the tangent to the breast curve which passes through the 
pixel. The tangent to the bottom half of the breast which passes through the point (Zq, 0) that 
is required passes through the point: 
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and thus: 



^tangent ^2 j, 



z 



V ^ ^ingent ) 

So, that for 0 < 5 < ffu,^,:, : h/a/dj = 0. Case B2 in Figure 23 reveals that for some of the 
5 angles of interest the x-ray path is pureiy through the breast edge. Now the two intersection 
points of the ray with the breast edge semi-circle are needed, i.e.(Zin,encci i, ^kamct i) and 

The two Z coordinates are given bv: 



m{a - c)± ^[{fTra^ - -f lac) 



'intersect 1 -f |77" 

10 with m -tan 0, a -H/2 and c = -tau 9Z,j, And thus: 



^^fai ) ~ ^intersect I '^imcrseci 2 ) ( ^ntereect 1 intersect 2 ) 

Using these thicknesses me basis extra-focal equation (7) can be rewritten as: 



0 



15 ] T(^)d^+ Je''*^''^'"^-'''''Y(^)^^ 



(8) 



where is the angle of the tangent to the breast edge curve and for cases Al and A2 is 
zero. 

In order to finaily compute the extra-focal radiation from equation (S) it is necessary 
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to estimate the constant The location of the projected breast edge is assumed and 

since the breast is assumed to be all fat the location of a curve which has /?,-„,=0 and hf^r^H 
(the inner curve in Figure 3). Thus the primary and scatter components along the curve can 
be estimated and by 

5 subtracting those values away from the total energy the extra-focal component along the 
curve can be found. Knowing just one value for the extra-focal component allows 
computation of the constant B"'^ extra as follows: 

Let the curve wich =0 and hj,~H be C/„„,=a. Since the location of that curve is 
known from the breast edge an estimate of the extra-focal component there is: 

10 V (;c,z) 6 Q,=o • Er''ix,z) = E';i{x,z) - E^M - E';\x,z) 

The scatter esumats comes from step (7) of the method, comes from the image 

and comes from knowing that /i,„=0 and /7/^,=H. The average of £?"^ along the C«„^o 
curve and Equation (S) v/iih 9 u„,i, = 0 to determine the constant: 



''extra /T/2 



0 



15 Thus the extra-focal component and the percentage of the total radiation that is 

exura-focal can be computed. Figure 24 shows two examples of the extra-focal component 
computed in this way. They show that, as expected, the extra-focal radiation is high at the 
breast edge but falls rapidly inside the breast falling to zero if the breast is large enough. 
Incorporating all the assumptions above in examples, the percentages computed 

20 indicate that extra-focal radiation is of ti-e order of 6-10% of the total radiation. This is on the 
low side of the estimate made in other ways which state that extra-focal component can make 
up to 15% of the total radiation, but is thought to be reasonable. 

Calculation of 

25 



The hi„, represe.itation can now be found by applying steps (9) and (10) to produce 
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surfaces such as that shown in Figure 5. In practice, steps (6) to (10) are used to calculate a 
first value of and then these steps are repeated with improved assumptions about the 
breast thickness H untii the best value of /?y„, is found. 
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CLAEVIS 

L A method of correcting a digital representation of an x-ray image for 
5 degradation caused by a digitiser, the digital representation being a measurement of the image 
density for each of a plurality of pixels of the image, the measurement being obtained by 
illuminating the image, measuring the attenuation of the light by the image and calculating 
from the attenuated light values the image density, the method comprising the steps of: 

calculating from the pixel values in the digital representation the value of the intensity 
10 of said attenuated light corresponding to each pixel; 

applying to the Intensity values a modulation transfer function to correct for the 
degradation introduced by the digitiser; and 

converting the corrected intensity values back into values representative of the image 
intensity. 

15 

2. A method according to claim 1, wherein the modulation transfer function 
corrects for blur and the digitising apenure. 

3. A method according to claim 1 or 2, wherein the intensity Ti of the attenuated 
20 light is calculated as: r=// x 10'^, where is the intensity of the illumination of the image and 

D is the value of the image density in the digital representation, and the corrected intensity 

values Ti' are converted back into values D' representative of the image 
intensity by: D' = log(I/Ti'). 

25 4. Apparatus for correcting a digital representation of an x-ray image for 

degradation caused by a digitiser, the digital representation being a measurement of the image 
density for each of a plurality of pixels of the image, the measurement being obtained by 
illuminating the image, measuring the a::enuaticn of the light by the image and calculating 
from the attenuated lig.:L values the image density, the apparatus comprising: 

30 calculating means for calculating from the pixel values in the digital representation 

the value of the intensity of said attenuated light corresponding to each pixel; 
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correction means for applying to the intensity values a modulation transfer function to 
correct for the degradation introduced by the digitiser; and 

conversion means for converting the corrected intensity values back into values 
representative of the image intensity. 

5 

5. Apparatus according to claim 4, wherein the conversion means corrects for 
blur and the digitising aperture. 

6. Apparatus according to claim 3 or 4, wherein the calculation means calculates 
10 intensity T/ of the attenuated light as: T:^Ii x 10'^, where I, is the intensity of the illumination 

of the image and D is the value of the image density in the digital representation, and the 
conversion means converts the corrected intensity values T/' back into values D' 
representative of the image intensity by: D ' = log(I/Ti 

15 7. A computer program storage medium readable by a computer system and 

encoding a computer program for controlling a computer to correct a digital representation of 
an x-ray image for degradation caused by a digitiser, the digital representation being a 
measurement of the im.age density for each of a plurality of pixels of the image, the 
measurement being obtained by illuminating the image, measuring the attenuation of the light 
20 by the image and calculating from the attenuated light values the image density, by a method 
comprising the steps of: 

calculating from the pixel values in the digital representation the value of the intensity 
of said attenuated light corresponding to each pixel; 

applying to the intensity values i modulation transfer function to correct for the 
25 degradation introduced by the digitiser; and 

convening the corrected intensify values back into values representative of the image 
intensity. 

8. A computer program storage medium according to claim 7, wherein the 
30 modulation transfer function corrects for blur and the digitising aperture. 
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9. A computer program storage medium according to claim 7 or 8, wherein the 
intensity Ti of the attenuated light is calculated as: Ti=I( x W^, where I, is the intensity of the 
illumination of the image and D is the value of the image density in the digital representation, 
and the corrected intensity values T/ are converted back into values D' representative of the 

5 image intensity by: D ' = log(I/T r 

10. In a method of x-ray imaging using an intensifying screen to receive x-rays 
and emit light to be recorded on an x-ray film, a method of calculating from the film image 
density the energy which was imparted to the intensifying screen, comprising the steps of: 

10 calibrating the film and intensifying screen by measuring the response of the film and 

intensifying screen to a plurality of different intensities of received x-rays; 

fitting a theoretical model of the expected response to the measured response, and 
using the fitted theoretical model to calculate the imparted energy from the image density. 

15 11. A method according to claim 10, wherein the theoretical model is a serpentine 

curve of the form: 

x^y^ a'y- b^x = 0, 

20 where a and b are constants, x is the logarithm of the energy imparted to the 

intensifying screen and y is the image density. ' ... ' 

12. A method according to claim 10 or 1 1, wherein the response of the film and 
intensifying screen to a plurality of different intensities of received x-rays is measured by 

25 exposing the film-screen combination to ::-rays through a lucite step wedge and measuring 

the image density procuced by the exposure through the different steps of the wedge. 

13. Apparatus for calculating from image densities of an x-ray film image the 
energy which was imparted to an intensifying screen used to receive x-rays and emit light to 

30 be recorded on the x-ray film, the apparatus comprising:- 

fitting means for fitting a theoretical model of the expected response of the x-ray film 
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and intensifying screen to a plurality of different intensities of received x-rays to a nieasured 
response, and 

calculation means for using the fitted theoretical model to calculate the imparted 
energy from the image density. 

5 

14. Apparatus according to ciaim 13, wherein the theoretical model is a serpentine 

curve of the form: 

x'y - a'y - b^x = 0, 

10 

where a and b arc: constants, x is - he logarithm of the energy imparted to the 
intensifying screen and y is the image density. 

15. A computer program storage m.edium readable by a computer system and 

1 5 encoding a computer program for controlling a computer to calculate from image densities of 
an x-ray film the energy which was imparted to an intensifying screen used to receive x-rays 
and emit light recorded on the x-ray fiin;, by a method comprising the steps of: 

fitting a theoretical model of the expected response of the x-ray film and intensifying 
screen to a plurality of different intensities of received x-rays to a measured response, and 

20 using the fitted theoretical model to calculate the imparted energy from the image density. 

16. A computer program storage medium according to claim 15, wherein the 
theoretical model is a serpentine curve of the form: 

25 x^y a'y - b'x - 0, 

where a and h e constants, x is :he logarithm of the energy imparted to the 
intensifying screen and y is the image c--jnsity. 

30 17. In a mernod of x-ray imaging in which an intensifying screen is used to 
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receive x-rays and emit light to be recorded on an x-ray, the image recorded on the film is 
digitised to give a pixelised digital representation of the image density on the film, a method 
of enhancing the digital representation of the image to remove the contribution thereto of 
glare from the intensifying screen comprising the steps of: 
5 converting the digital representation of image density into a representation of the 

energy impaned to the screen; and 

deconvolving the representation of the energy imparted to the screen using a 
weighting mask defining the point spread function for the intensifying screen to remove the 
contribution thereto of glare from the intensifying screen. 

10 

18. A method according to claim 1 7, wherein the weighting mask defines the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels weighted by distance from the 
absorption site to the pixel. 

15 

19. A method according to ciaim 18 , wherein the weighting mask defines the 
contribution to the energy imparted to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels further weighted by die depth of 
the absorption site in the intensifying screen. 

20 

20. A method of detecting fihm-screen shot noise in a x-ray image in which an 
intensifying screen is used to receive x-rays and to emit light to be recorded as an image on a 
x-ray film, the image recorded on the film being digitised to give a pixelised digital 
representation thereof . /hich is then converted into a representation of the energy imparted to 

25 the screen, the method comprising removing from that representation glare from the 

inter.sifying screen in accordance with crie method of claim 17, and defining as film-screen 
shot noise regions of the image in whicii the resultant value for energy imparted to the screen 
is infeasible. 

30 21. A meth r d according to claim 20, wherein regions of the image in which 

resultant values for er. -rgy imparted to :he screen having negative values are defined as film- 
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22. A method according to claim 20 or 21, wherein the image is a mammogram. 

5 23 . Apparatus for enhancing the digital representation of an x-ray image obtained 

by using an intensifying screen to receive x-rays and emit light to be recorded on an x-ray 
film and digitising the film to give a pixelised digital representation of the image density 
thereon the apparatus comprising: 

conversion means for convening: the digital representation of image density into a 
10 representation of the energy imparted to [he screen; and 

deconvolution means for decon volving the representation of the energy imparted to 
the screen using a weighting mask defining the point spread function for the intensifying 
screen to remove the contribution thereto of glare from the intensifying screen. 

15 24. Apparatus according to claim 23, wherein the weighting mask defines the 

contribution to the energy impaned to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels weighted by distance from the 
absorption site to the pixel. 

20 25. Apparatus according to claim 24 , wherein the weighting mask defines the 

contribution to the energy impaned to the screen caused by the emission of light from the 
point of absorption of an x-ray towards neighbouring pixels further weighted by the depth of 
the absorption site in the intensifying screen. 

25 26. Apparatus for detecting film-screen shot noise in a x-ray image obtained by 

using an intensifying screen to receive x-rays and to emit light to be recorded as an image on 
an x-ray film, the image recorded on the film being digitised to give a pixelised digital 
representation thereof which is then ccriverted into a representation of the energy imparted to 
the screen, the apparatus comprising apparatus according to claim 23, 24 or 25 for removing 

30 frorr. that representation glare from the intensifying screen, and means for defining as film- 
screen shot noise regions of the image in which the resultant value for energy imparted to the 
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27. Appararus according to claim 26, wherein regions of the image in which 
resultant values for energy imparted to the screen having negative values are defined as film- 

5 screen shot noise. 

28. Apparatus according to claim 26 or 27, wherein the image is a mammogram, 

29. A computer program storage medium readable by a computer system and 

10 encoding a computer program for controlling a computer to enhance the digital representation 
of an x-ray image obtained by using an intensifying screen to receive x-rays and emit light to 
be recorded on an x-ray film, the image recorded on the film being digitised to give a 
pixelised digital representation of the image density thereon, the removal from the digital 
representation of the image of the contribution thereto of glare from the intensifying screen 

1 5 comprising the steps of: 

converting the digital representation of image density into a representation of the 
energy impaned to the screen; and 

deconvolving the representation of the energy impaned to the screen using a 
weighting mask defining the point spread function for the intensifying screen to remove the 

20 contribution thereto of glare from the intensifying screen. 

30. A computer program storage medium according to claim 29, wherein the . 
weighting mask defines the contribution to the energy imparted to the screen caused by the 
emission of light from the point of absorption of an x-ray towards neighbouring pixels 

25 weighted by distance from the absorption site to the pixel. 

31. A computer program storage medium according to claim 30 , wherein the 
weighting mask defines the contributior. lo the energy imparted to the screen caused by the 
emission of light from the point of absorption of an .x-ray towards neighbouring pixels further 

30 weighted by the depih of the absorption site in the intensifying screen. 



10 
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32. A computer program storage medium readable by a computer system and 
eiicoding a computer program for contrciiing a computer to detect film-screen shot noise in a 
x-ray image obtained by using an intensifying screen to receive x-rays and to emit light to be 
recorded as an image on an x-ray film, the image recorded on the film being digitised to give 
a pixelised digital representation thereof which is then converted into a representation of the 
energy imparted to the screen, the detection of film-screen shot noise comprising removing 
from that representation glare from the intensifying screen in accordance with claim 29, and 
defining as film-screen shot noise regions of the image in which the resultant value for 
energy imparted to the screen is infeasib-e. 



33. A computer program storage medium according to claim 32, wherein regions 
of the image in which resultant values f:. : energy impaned to the screen having negative 
values are defined as film-screen shot noise. 

15 34. A computer program storage medium according to claim 32 or 33, wherein the 

image is a mammogram. 

35. A method of calculating from a mammogram the compressed thickness of the 
imaged breast, comprising the step of delimiting in the mammogram the region 

20 corresponding to the part of the breast which is compressed from the region corresponding to 
the uncompressed breast edge by detecnng the smoothness of curv es of equal intensity in the 
mammogram. . > 

36. A method according to Ciaim 35, further comprising the steps of detecting the 
25 smooth equal intensity curve which is positioned furthest from the breast edge, calculating 

from the intensities along that curve the thickness of fat in the breast which gives those 
intensities, and calcuiating the thicknesi; of the compressed breast from that thickness of fat. 



30 



37. A mechod according to claim 35 or 36, comprising as an initial step the 
conversion of the mammogram into the representation representing the thicknesses of fat 
and interesting tissue in regions of the breast contributing to each point in the mammogram. 
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38. A merhod according to claim 37, wherein the conversion is performed on the 
basis of an underestimate of the compressed breast thickness, the conversion being iterated 
with successively changing estimates of :he comxpressed breast thickness until a smooth 
delimitation between the region correspcnding to the part of the breast which is compressed 

5 and the region corresponding to the uncompressed breast edge is found. 

39. A method according to any one of claims 35 to 38, wherein a predefined 
smoothness threshold is set to detect the smooth curves. 

10 40. A method of calculating cne contribution to a mammogram of extra-focal 

radiation by calculating the compressed :::ickne5S of the imaged breast according to the 
method of any one of claims 35 to 39, cal culating from the x-ray exposure and an estimate of 
the scattering of radiation the expected ir.iensicy in the mammogram along the delimitation 
between the region corresponding to the paa of the breast which is compressed and the region 

15 corresponding to the uncompressed breas: edge, and comparing the expected intensity to the 
actual intensity in the mammogram. 

41. Apparatus for calculating from a mammogram the compressed thickness of 
the imaged breast, comprising means for delimiung in the mammogram the region 

20 corresponding to the part of the breast vvliich is compressed from the region corresponding to 
the uncompressed breast edge by detecting the smoothness of curves of equal intensity in the 
mammogram. 

42. Apparatus according to claim 41, funher comprising detection means to detect 
25 the smooth equal intensity curve which in positioned furthest from the breast edge, and 

calculating means fcr calculating from intensities along that curve the thickness of fat in 
the breast which gives those intensities^ ^md for calculating the thickness of the compressed 
breast from that thic>:ness of fat. 

30 43. Apparatus according to :':.i.m 41 or 42, comprising conversion means for 

converiing the mammogram into the /?;,,. representation representing the thicknesses of fat and 
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contributing to each point in the mammogram. 



44. Apparatus according to ciaim 43, wherein the conversion means first 
underestimates the compressed breast thicicness, and iterates the conversion with successively 

5 changing estimates of the compressed breast thickness until a smooth delimitation between 
the region corresponding to the part of the breast which is compressed and the region 
corresponding to the uncompressed breas: edge is found. 

45. Apparatus according to aay one of claims 41 to 44, wherein a predefined 
10 smoothness threshold is set to detect the smooth curves. 

46. Apparatus for calculating the contribution to a mammogram of extra-focal 
radiation comprising apparatus according to any one of claims 41 to 45 for calculating the 
compressed thickness of the imaged breas r, and further comprising expected intensity 

15 calculating means fcr calculating from ta e x-ray exposure and an estimate of the scattering of 
radiation the expected intensity in the mammogram along the delimitation between the region 
corresponding to the part of the breast which is compressed and the region corresponding to 
the uncompressed breast edge, and comparing means for comparing the expected intensity to 
the actual intensity in the mammogram. 

20 

47. A computer program storage medium readable by a computer system and 
encoding a computer program for contro- ling a computer to calculate from a mammogram the 
compressed thickness of the imaged breast by a method comprising the step of delimiting in 
the mammogram the region corresponding to the part of the breast which is compressed from 

25 the region correspoi\ding to the uncompressed breast edge by detecting the smoothness of 
curves of equal intensity in the mammci^ram. 

43. A computer program stcr-ge medium according to claim 47, wherein the 
method further comprises the steps of ce:acting the smooth equal intensity curve which is 
30 positioiied funhest from the breast edge, calculating from the intensities along that curve the 
thickness of fat in t-ie breast which gives :hose intensities, and calculating the thickness of the 
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compressed breast from that tfaicloiess of fat. 

49. A computer program storage medium according to claim 47 or 48, wherein the 
method fjriher comprises as an initial step the conversion of the mammogram into the A,-^ 

5 representation representing the thicknesses of fat and interesting tissue in regions of the breast 
contributing to each point in the mammogram. 

50. A coinputer program store.ge medium according to claim 49, wherein the 
conversion is performed on the basis ofzn underestimate of the compressed breast thickness, 

10 the conversion being iterated with successively changing estimates of the compressed breast 
thickness until a smooth delimitation betv/een the region corresponding to the part of the 
breast which is compressed and the regio:i corresponding to the uncompressed breast edge is 
found. 

15 5 : ^ A computer program storage medium according to any one of claims 47 to 50, 

wherein a predefined smoothness threshold is set to detect the smooth curves. 

52. A computer program storage medium readable by a computer system and 
encoding a computer program for controlling a computer to calculate the contribution to a 

20 mammogram of extra-focal radiation by calculating the compressed thickness of the imaged 
breast according to ihe method of any one of claims 35 to 40, calculating from the x-ray 
exposure and an estimate of the scattering of radiation the expected intensity in the 
mammogram along the delimitation between the region corresponding to the part of the 
breast which is compressed and the regici: corresponding to the uncompressed breast edge, 

25 and comparing the expected intensity to tae actual intensity in the mammogram. 

53 . A method of processing a inammogram to produce an hi„f representation 
represe:- :;ag the thicknesses of fat and irieresting tissue in regions of the breast contributing 
to each coint in the mammogram, the mammogram being produced by using an intensifying 

30 screen to receive x-rays transmitted thrcgh the breast and to emit light to be recorded as the 
mammogram on an x-ray film, the mam:::ogram being digitised by a digitiser to obtain a 
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digital representation thereof comprising p ixel values representing the image intensities in the 
mammogram, the method comprising the steps of: 

correcting the digital representation for degradation caused by the digitiser according 
to the method of any one of claims I to 3 , 

converting the values in the digita. representation to values representing the energy 
which was imparted to the intensifying screen according to the method of any one of claims 
10 to 12; 

compensating the converted values for glare from the intensifying screen according to 

the method of any one of claims 12 to 22; 

funher compensating the values for variations in intensity of the x-ray beam incident 

on the breast; 

calculating the compressed thickness of the imaged breast according to the method of 

any one of claims 35 to 39; 

calculating the contribution to the mammogram of extra-focal radiation according to 

die method of claim 40; 

s'jbcracting from the compensated values the extra-focal radiation and an estimate of 
the scattered radiation to give enhanced measured values of the energy which was imparted to 
the intensifying screen; and 

calculating the A,„, representation by equating the enhanced measured values with 

predicted values. 

54. A method according to claim 53 , wherein the predicted values are calculated 
from the intensity of x-rays incident on the breast and the absorption characteristics of 
physical elements of die imaging system. 

: 5. Apparatus for processing a mammogram to produce an /7,„, representation 
represen ting the thicknesses of fat and in teresting tissue in regions of the breast contributing 
to each -oint in the mammogram, the m-t.nmogram being produced by using an intensifying 
screen i ; receive x-rays transmitted thro-.gh the breast and to emit light to be recorded as the 
mammc gram on an x-ray film, the mamiviogram being digitised by a digitiser to obtain a 
dieital r . nresentation thereof comprising pixel values representing the image intensities in the 
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mammogr-m, the apparatus comprising: 

apparatus for correcting the digital representation for degradation caused by the 
digitiser according to any one of claims 4 to 6; 

apnaratus for calculating from the values in the digital representation values 
representir.g the energy which was imparted to the intensifying screen according to claim 13 
or 14; 

apparatus for compensating the converted values for glare from the intensifying 
screen according to any one of claims 23 to 28; 

means for further compensating the values for variations in intensity of the x-ray 

beam inc: .lent on the breast; 

ap::aratus for calculating the compressed thickness of the imaged breast according to 

any one of claims 41 to 45 

ap :;aratus for calculating the contribution to the mammogram of extra-focal radiation 

according to claim 46; 

subtraction means for subtracting from the compensated values the extra-focal 
radiation i.nd an estimate of the scattered radiation to give enhanced measured values of the 
energy vviuch was imparted to the intensifying screen; and 

hi,, calculating means for calculating the representation by equating the enhanced 
measurec values with predicted values. 

56 . A computer program storage medium readable by a computer system and 
encoding i computer program for controlling a computer to process a mammogram ta 
produce a:; representation representing the thicknesses of fat and interesting tissue in 
regions c ': the breast contributing to each point in the mammogram, the mammogram bemg 
producer using an intensifying screen ■ o receive x-rays transmitted through the breast and 
to emit li :ht to be recorded as the mammogram on an x-ray film, the mammogram being 
digitised -y a digitiser to obtain a digital representation thereof comprising pixel values 
represen: the image intensities in the ir.ammogram, by a method comprising the steps of: 

cc-recting the digital representaticf: for degradation caused by the digitiser according 
to the me nod of any one of claims 1 to 3; 

c -liverting the values in the digita; representation to values representing the energy 
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which was imparted to the intensifying screen according to the method of any one of claims 
10 to 12; 

compensating the converted values for glare from the intensifying screen according to 
the method of any one of claims 12 to 22; 
5 further compensating the values for variations in intensity of the x-ray beam incident 

on the breast; 

calculating the compressed thickness of the imaged breast according to the method of 

any one of claims 35 to 39; 

c:.lculating the contribution to the mammogram of extra-focal radiation according to 

10 the metivod of claim 40; 

subtracting from the compensated values the extra-focal radiation and an estimate of 
the scattered radiation to give enhanced n:easured values of the energy which was imparted to 
the intensifying screen; and 

calculating the /7,„, representation by equating the enhanced measured values with 

15 predicted values. 

57. A method of detecting microcalcifications in a breast from a mammogram of 
the breast, comprising the steps of processing the mammogram to produce an 
represeniiiion according to the method of claim 53 or 54, converting the value of /2,-„, for a 
20 candidate region into a value representative of the volume of interesting tissue and 

thresho!::ing the volume values to detect as microcalcifications areas of the mammogram 
where the volume value exceeds a threshold. 

52. A method according to claim 57 wherein the value of the volume of interesting 
25 tissue is compared to an estimate of the actual volume of the candidate region. 

5 5. A method according to claim 57 or 58 wherein the volume values are 
calculated after subtraction from the h:„, value of an estimate of the background 

30 6-. . Apparatus for detecting microcalcifications in a breast from a mammogram of 

the brea'; , comprising apparatus for processing the mammogram to produce an A/^, 
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representa:ion according to claim 55, means for convening the value ofhi„, for a candidate 
region into a value representative of the volume of interesting tissue and means for 
thresholding the volume values to detect as microcalcifications areas of the mammogram 
where the volume value exceeds a threshold. 

5 

61 . Apparatus according to claim 60 further comprising means for comparing the 
value of tl;e volume of interesting tissue to an estimate of the actual volume of the candidate 
region. 

10 62. Apparatus according to claim 60 or 61 further comprising means for 

subtracting from the value an estimate of the background before calculating the 
volume values. 

63. A computer program storage medium readable by a computer system and 
15 encoding a computer program for controlling a computer to detect microcalcifications in a 

breast frcni a mammogram of the breast by a method comprising the steps of processing the 
raammogi am to produce an /7,„^ representation according to the method of claim 53 or 54, 
convening the value of for a candidate region into a value representative of the volume of 
interesting tissue and thresholding the volume values to detect as microcalcifications areas of 
20 the mammogram where the volume value exceeds a threshold. 

64. A computer program storage medium according to claim 63 wherein the value 
of the vc/ume of interesting tissue is compared to an estimate of the actual volume of the 
candidate region. 

25 

6i . A program storage medium according to claim 63 or 64 wherein the volume 
values ar j calculated after subtraction from the h:„, value of an estimate of the background 



30 



c . A method of processing a mammogram substantially as hereinbefore described 
with reference to and as illustrated in the accompanying drawings. 
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